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Abstract 

We reanalyze deep inelastic scattering data of BCDMS Collaboration by in- 
cluding proper cuts of ranges with large systematic errors. We perform also the 
CsJ ' fits of high statistic deep inelastic scattering data of BCDMS, ShAC, NM and BFP 

Collaborations taking the data separately and in combined way and find good agree- 
\ ment between these analyses. We extract the values of the QCD coupling constant 

, i ; u s {Mz) up to NLO level. 

; 1 Introduction 

The deep inelastic scattering (DIS) leptons on hadrons is the basical process to study 
\ the values of the parton distribution functions (PDF) which are universal (after choosing 
of factorization and renormalization schemes) and can be used in other processes. The 
accuracy of the present data for deep inelastic structure functions (SF) reached the level 
at which the Q 2 -dependence of logarithmic QCD-motivated terms and power-like ones 
may be studied separately (for a review, see the recent papers [p]] and references therein). 
In the present letter we sketch the results of our analysis |2| at the next-to-leading 
p order (NLO) of perturbative QCD for the most known DIS SF F 2 (x,Q 2 ) [] taking into 
account experimental data [g-0 of SLAC, NM, BCDMS and BFP Collaborations. We 
stress the power-like effects, so-called twist-4 (i.e. ~ 1/Q 2 ) contributions. To our purposes 
we represent the SF F 2 (x, Q 2 ) as the contribution of the leading twist part F 2 D {x, Q 2 ) 
described by perturbative QCD, when the target mass corrections are taken into account 
(and coincides with F 2 tw2 (x, Q 2 ) when the target mass corrections are withdrawn), and 
^ ■ the nonperturbative part ( "dynamical" twist-four terms) : 

F 2 (x,Q 2 ) = Ft u (x,Q 2 ) = Ff CD {x,Q 2 ) (l + **M), (1) 

where h^x) is magnitude of twist-four terms. 

Contrary to standard fits (see, for example, |§- [[K]]) when the direct numerical calcu- 
lations based on Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation [JTTJ] are 
used to evaluate structure functions, we use the exact solution of DGLAP equation for 
the Mellin moments ilC^Q 2 ) of SF F 2 tw2 (x,Q 2 ): 

M*(Q 2 ) = J x n ~ 2 F^(x,Q 2 )dx (k = full,pQCD,tw2, ...) (2) 

and the subsequent reproduction of F%(x, Q 2 ) at every needed Q 2 -value with help of 
the Jacobi Polynomial expansion method [|l|, |l3j (see similar analyses at the NLO level 
JI3|. n] and at the next-next-to- leading order (NNLO) level and above flTo| . 

In this letter we do not present exact formulae of Q 2 -dependence of SF F 2 which are 
given in 0. We note only that the moments M^ 2 {Q 2 ) at some Ql is theoretical input of 

^ere Q 2 = -q 2 and x = Q 2 /(2pq) are standard DIS variables, where q and p are photon and hadron 
momentums, respectively. 
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our analysis and the twist-four term h±{x) is considered as a set of free parameters (one 
constant h^Xi) per Xj-bin): h{ ree (x) = Yl!i=\ ^4(^1) > where / is the number of bins. 

2 Fits of F 2 : procedure 

Having the QCD expressions for the Mellin moments we can reconstruct the SF F^{x) 
as 

Afi max Tl 

F*' N ™*(x,Q 2 ) = x a (l-x) b £ 9f (*) J>f («,/3)MjV 2 (Q 2 ) , (3) 

n=0 j=0 

where 0°' 6 are the Jacobi polynomials Q and a, 6 are fitted parameters. 

First of all, we choose the cut Q 2 > 1 GeV 2 in all our studies. For Q 2 < 1 GeV 2 , the 
applicability of twist expansion is very questionable. Secondly, we choose quite large values 
of the normalization point Q^: our perturbative formulae should be applicable at the value 
of Ql- Moreover, the higher order corrections ~ a^(Ql) and ~ (a s (Q 2 ) — a s (Ql)) k (k > 2) 
should be less important at higher Ql values. 

We use MINUET program p[ for minimization of x 2 (^ 2 ) = \{F™ P - F* eor ) / AF* xp \ 2 . 



We consider free normalizations of data for different experiments. For the reference, we 
use the most stable deuterium BCDMS data at the value of energy E = 200 GeV (E 
is the initial energy lepton beam). Using other types of data as reference gives negligible 
changes in our results. The usage of fixed normalization for all data leads to fits with a 
bit worser \ 2 - 

3 Results of fits 

Hereafter we choose Qq = 90 GeV 2 (Qq = 20 GeV 2 ) for the nonsinglet (combine nonsinglet 
and singlet) evolution, that is in good agreement with above conditions. We use also 

3.1 BCDMS 12 C + H 2 + D 2 data 

We start our analysis with the most precise experimental data obtained by BCDMS 
muon scattering experiment at the high Q 2 values. The full set of data is 762 (607) points 
(for the bounded x range: x > 0.25). 

It is well known that the original analyses given by BCDMS Collaboration itself (see 
also Ref. 0) lead to quite small values a s (M|) = 0.113. Although in some recent papers 
(see, for example, || 20||) more higher values of the coupling constant a s (M|) have been 



observed, we think that an additional reanalysis of BCDMS data should be very useful. 
Based on study J21| we proposed in that the reason for small values of a s (M|) 



coming from BCDMS data was the existence of the subset of the data having large sys- 
tematic errors. We studied this subject by introducing several so-called V-cutsQ (see [0). 
Excluding this set of data with large systematic errors leads to essentially larger values of 
a s (M|) and very slow dependence of the values on the concrete choice of the Y-cut (see 
below). 

We use the following x-dependent Y-cuts: 

2 We note here that there is similar metho d JL6| ], based on Bernstein polynomials. The method has 
been used in the analyses at the NLO level in Jl7[ and at the NNLO level in 

3 Hereafter we use the kinematical variable Y = (E — E)/E , where E is scattering energies of lepton. 
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Figure 1: The study of systematics at 
different Y cut values in the fits based on 
nonsinglet evolution. The QCD analysis of 
BCDMS 12 C, H 2 , D 2 data (nonsinglet case) 
is given at x cut = 0.25 and Ql = 90 GeV 2 . 
The inner (outer) error-bars show statisti- 
cal (systematic) errors. 



Figure 2: The study of systematics at dif- 
ferent Y C ut values in the fits based on com- 
bine singlet and nonsinglet evolution. All 
other notes are as in Fig. 1 with two ex- 
ceptions: no a x cut and Ql = 20 GeV 2 . 
Moreover, the points Ny C ut = 1,2,3,4,5 
correspond the values N = 1,2,4,5,6 in 
the Table 1. 



when 0.3 < x < 0.4, y > 0.16 when 0.4 < x < 0.5 
when 0.5 < x < 0.6, y > Y cu u when 0.6 < x < 0.7 
when 0.7 < x < 0.8 (4) 
for the cuts at 0.5 < x < 0.8: 
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0.16 
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0.20 


0.20 


0.22 


0.22 


0.24 


0.25 



Table 1: The values of Y cut3} Y cut4 and Y cut5 . 



The systematic errors for BCDMS data were given || as multiplicative factors to be 
applied to F 2 (x,Q 2 ): f r , fb, fs, fd and fh are the uncertainties due to spectrometer res- 
olution, beam momentum, calibration, spectrometer magnetic field calibration, detector 
inefficiencies and energy normalization, respectively. For this study each experimental 
point of the undistorted set was multiplied by a factor characterizing a given type of un- 
certainties and a new (distorted) data set was fitted again in agreement with our procedure 
considered in the previous section. The factors (f r , fb, f s , fd, fh) were taken from papers [|] 
(see CERN preprint versions in ||). The a s values for the distorted and undistorted sets 
of data are given in the Figs. 1 and 2 (for the cases of nonsinglet and complete evolutions, 
respectively) together with the total systematic error estimated in quadratures. 

From the Figs. 1 and 2 we can see that the a s values are obtained for iV = 1 4- 6 of 
Y cu t3, Y cut 4 and Y cut5 are very stable and statistically consistent. The case iV = 6 of the 



y > 0.14 

y > Y cut3 

y > Ycutb 

and several iV sets 
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Table 1 reduces the systematic error in a s by factor 1.8 and increases the value of a s , 
while increasing the statistical error on the 30%. 

After the cuts have been implemented (we use the set N = 6 of the Table 1), we have 
590 (452) points (for the bounded x range: x > 0.25). Fitting them in agreement with 
the same procedure considered in the previous Section, we obtain the following results: 

from fits, based on nonsinglet evolution (i.e. when x > 0.25): 

a s (M 2 z ) = 0.1153 ± 0.0013 (stat) ± 0.0022 (syst) ± 0.0012 (norm), 

from fits, based on combined singlet and nonsinglet evolution: 

a s (M 2 z ) = 0.1175 ± 0.0014 (stat) ± 0.0020 (syst) ± 0.0011 (norm), (5) 

where hereafter the symbol "norm" marks the error of normalization of experimental data. 

The results are agree each other within considered errors. In Ref. H we have also 
analyzed the combine SLAC, NM and BFP data and found good agreement with (|5]). 
So, we have a possibility to fit together all the data. It is the subject of the following 
subsection. 

3.2 SLAC, BCDMS, NM and BFP data 

After these V-cuts have been incorporated (with N = 6) for BCDMS data, the full set of 
combine data is 1309 (797) points (for the bounded x range: x > 0.25). 

To verify the range of applicability of perturbative QCD, we analyze firstly the data 
without a contribution of twist-four terms, i.e. when F 2 = F 2 . We do several fits 
using the cut Q 2 > Q 2 cut and increase the value Q 2 cut step by step. We observe good 
agreement of the fits with the data when Q 2 cut > 10 4- 15 GeV 2 (see the Figs. 3 and 4). 
Later we add the twist-four corrections and fit the data with the standard cut Q 2 > 1 
GeV 2 . We have find very good agreement with the data. Moreover the predictions for 
a s (M z ) in both above procedures are very similar (see the Figs. 3 and 4). The results of 
the fits are compiled in Summary (see Eqs. (§|)-(0)). 

4 Summary 

We have demonstrated several steps of our study of the Q 2 -evolution of DIS structure 
function F 2 fitting all modern fixed target experimental data. 

From the fits we have obtained the value of the normalization a s (M z ) of QCD coupling 
constant. First of all, we have reanalyzed the BCDMS data cutting the range with large 
systematic errors. As it is possible to see in the Fig. 1, the value of a s (M z ) rises strongly 
when the cuts of systematics were incorporated. In another side, the value of a s (M z ) 
does not dependent on the concrete type of the cut within modern statistical errors. 

We have found that at Q 2 > 10-rl5 GeV 2 the formulae of pure perturbative QCD (i.e. 
twist-two approximation together with target mass corrections) are in good agreement 
with all data. [] The results for a s (M z ) are very similar (see ||) for the both types of 

4 We note that at small x values, the perturbative QCD works well starting with Q 2 = 1.5-^2 GeV 2 and 
higher twist corrections are important only at very low Q 2 : Q 2 ~ 0.5 GeV 2 (see |2^| and references 
therein). As it is was observed in |24|, ^5) (see also discussions in |2^, ^6|) the good agreement between 
perturbative QCD and experiment seems connect with large effective argument of coupling constant at 
low x range. 
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Figure 3: The values of a s {M\) and x 2 
at different Q 2 -values of data cuts in the 
fits based on nonsinglet evolution. The 
black (white) points show the analyses of 
data without (with) twist-four contribu- 
tions. Only statistical errors are shown. 



Figure 4: The values of a s (M|) and x 2 
at different Q 2 -values of data cutes in the 
fits based on combine singlet and nonsin- 
glet evolution. All other notes are as in 
Fig. 3. 



analyses: ones, based on nonsinglet evolution, and ones, based on combined singlet and 
nonsinglet evolution. They have the following form: 

• from fits, based on nonsinglet evolution: 

a s (M 2 z ) = 0.1170 ± 0.0009 (stat) ± 0.0019 (syst) ± 0.0010 (norm), (6) 

• from fits, based on combined singlet and nonsinglet evolution: 

a s (M 2 z ) = 0.1180 ± 0.0013 (stat) ± 0.0021 (syst) ± 0.0009 (norm), (7) 

When we have added twist-four corrections, we have very good agreement between 
QCD (i.e. first two coefficients of Wilson expansion) and data starting already with 
Q 2 = 1 GeV 2 , where the Wilson expansion should begin to be applicable. The results for 
a s (M|) coincide for the both types of analyses: ones, based on nonsinglet evolution, and 
ones, based on combined singlet and nonsinglet evolution. They have the following form: 

• from fits, based on nonsinglet evolution: 

a s (M 2 z ) = 0.1174 ± 0.0007 (stat) ± 0.0019 (syst) ± 0.0010 (norm), (8) 



5 



• from fits, based on combined singlet and nonsinglet evolution: 



a s (M|) = 0.1177 ± 0.0007 (stat) ± 0.0021 (syst) ± 0.0009 (norm), (9) 

Thus, there is very good agreement (see Eqs. (||), (|7|), (||) and @) between results 
based on pure perturbative QCD at quite large Q 2 values (i.e. at Q 2 > 10-r- 15 GeV 2 ) and 
the results based on first two twist terms of Wilson expansion (at Q 2 > 1 GeV 2 , where 
the Wilson expansion should be applicable). 

We would like to note that we have good agreement also with the analysis J20| of 
combined HI and BCDMS data, which has been given by HI Collaboration very recently. 
Our results for a s (M|) are in good agreement also with the average value for coupling 



constant, presented in the recent studies (see || £7], [18], |28fl and references therein) and 
in famous Altarelli and Bethke reviews [p9|j . 

The last result (||) based on all data with Q 2 > 1 GeV 2 can be considered as "best 
value" for the coupling constant a s {M^) coming in our analysis. 
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